On solutions of the nonequilibrium x-ray edge problem 
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We rediscuss a nonequilibrium x-ray edge problem which in recent publications led to discrepancies 
between the results of the perturbative and of an extended Nozieres-De Dominicis approach. We 
show that the this problem results from an uncritical separation of momenta of the scattering 
potential, and we propose a corrected Nozieres-De Dominicis solution. 



I. INTRODUCTION 

We address a question raised recently by Combescot 
and Rouletn (CR) in the context of the nonequilibrium 
x-ray edge, problem. The latter was studied some years 
ago by Ngp who developed an extension of the well-known 
Nozieres and De DominicisB (ND) technique. In the late 
60s, ND devised this method for the computation of the 
x-ray absorption spectrum. They provided with it an el- 
egant description of the shape of the edge singularity in 
terms of the scattering phase shift at the Fermi surface. 
The approach is based on some subtle approximations, 
though, and thorough check of the condition of validity 
is required. Such a check is provided, for instance, by 
comparison with the perturbative pjaxquet diagram re- 
sult of Roulet, Gavoret and NozieresEfB, and the validity 
of ND's solution for the original x-ray problem can in this 
sense be proven. For the nonequilibrium case, however, 
CR revealed an inconsistency by comparing a similar par- 
quet diagram summation to Ng's result, yet they did not 
provide an explanation of its origin. In the present paper, 
we show that the discrepancies arise from an uncritical 
use of the decoupling of the momenta of the scattering 
potential. Further, we show how ND's method can be 
modified to yield a result coinciding with the perturba- 
tive one. This is illustrated on CR's academic problem 
for which it turns out that Ng's multi-channel extensions 
are not requiredB. 



II. NONEQUILIBRIUM X-RAY PROBLEM 

The nonequilibrium system to be discussed consists of 
two Fermi seas n = 1,2, referred to as subbands which are 
characterized by the fixed chemical potentials fi± < \i2 
and the energies £1,62 at the bottom of the subbands. 
We assume that [i\ > £2, so there is a nonzero overlap 
between the bottom of the subband 2 and the Fermi sur- 
face Hi. A core state with energy Ed and an infinite mass 
is assumed to exist below the subbands. The qualitative 
situation is shown in Fig. [l]. 

If both subbands were isolated, the x-ray absorption 
spectrum would be the result of a simple superposition 
of two independent x-ray problems, and two edge singu- 
larities corresponding to [i\ and (X2 could be observed. 
The mixture between the two subbands, however, allows 



for a new physical process: Particle-hole excitations with 
the particle in one subband and the hole in the other one 
become possible with energies that may be positive or 
negative (see Fig. |l|). The form of the absorption spec- 
trum is changed considerably. The edge singularities are 
broadened, the core-hole acquires a finite life-time, and 
absorption below the threshold becomes possible, where 
the missing energy is compensated by such particle-hole 
excitations of negative energy. 

The Hamiltonian of the system is assumed to be given 

by 
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kn 
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where n = 1,2 labels the subbands. Conduction elec- 
trons with momenta k in the subbands n are created and 
annihilated by ct and Ckn, while b^ and b operate on 
the core state. The interaction with the photon field is 
described by the operator 



(2) 



kn 



For simplicity, spin indices have been dropped. 



III. PERTURBATIVE APPROACH OF 
COMBESCOT AND ROULET 



The problem 
ans (|l|) and (f 



as it is defined by the Hamiltoni- 
above, depends on the parameters 
i 1 ^fefe 1 ' i ^4 2 i/~ an d Vkk' ■ With respect to the original x- 
ray problem the new physics must enter with the subband 
couplings Vuu- To focus on the new effects we consider 
with CR the academic problem of V k l y = V k 22 , = A| = 0. 
This means that a photon can excite the core electron 
only into the subband n — 1, and particle-hole excita- 
tions are due only to particles tunneling between both 
subbands, as indicated in Fig. |]. 

The perturbative approach of CR consists in assuming 
V 12 and in expanding the response function 
5(£!), which is the Fourier transform of 
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in powers of 
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The imaginary part of <S(0) yields the absorption rate 
A(Cl). CR performed the sum over the most singular 
parquet diagrams and found an absorption rate of the 
x-ray problem which is, close to the threshold energy 
^0 = Mi ~ £d, of the form 



A(n) « Im 



2g u \ 



2(/ 



(5) 



for £ a cutoff of the order of the band width, 9 the usual 
step function, and g the effective coupling constant 



9 = vnV(v 2 i - iirv 2 x). 



(6) 



Here we have used the notation v nn i = v n (/j,' n ) for the 
density of states in the subband n at the energy \i n i , and 
v nn i for the integral 
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with P denoting the principal value. The approximation 
(||) holds for small values of the coupling constant g. 

CR compared this result with Ng's extended ND ap- 
proach. Their criticism refers to the appearance of the 
density of states V22 = ^2(^2), and not, as in the per- 
turbative result above, of the density of states at 
^21 = ^2(^1)- The presence of z/22 implies a coupling of 
particles at the respective energies \i\ and /j,2, whereas 
the zero energy inter-subband fluctuations at \i\ are en- 
tirely absent. This is definitely unphysical as soon as 
A'2 — Mi becomes large. The origin of this discrepancy 
remained unresolved in CR's paper and is discussed in 
the next section. 



IV. ORIGIN OF THE UNPHYSICAL 
BEHAVIOR 

The problems with ND's approach arise from an un- 
critical use of the momentum decoupling 
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As CR have noticed with the perturbative results above, 
the important particle-hole excitations occur in the low 
energy sector with the particle and hole energies close 
to the Fermi surface Ml- In the single-band equilibrium 
x-ray problem, the restriction to this energy sector can 
be assured by the decoupling (^). As it turns out, the 
explicit form of the is of little importance because 
the form of the free propagator is controlled by the dis- 
continuity at the Fermi surface (see Ref. 3, Eq. (33) 
and below) which naturally imposes the low-energy re- 
strictions. There is consequently no risk in assuming the 
potential V^k' to be constant in the whole band as long 
as one introduces the adequate cutoff £0 imposed by the 
band width. 

In the noncquilibrium case a similar "naive" splitting 
of the form (0) leads to the unphysical results mentioned 



above. If one assumes the potential to be constant over 
the whole band widths, the free propagators which are 
integrated over the u k are controlled by the discontinu- 
ities at the Fermi surfaces of their respective subbands 
analogously to the original x-ray problem. The coupling 
of the subbands then leads to the mixture of the energy 
sectors close to the different Fermi surfaces, and there- 
fore to large energy particle-hole excitations of the order 
of fJ-2 — Mii whereas the inter-subband fluctuations close 
to Hi and /i 2 respectively are absent. From CR's pertur- 
bative approach, however, we know that the latter are 
the relevant processes and cannot be neglected. 

In the sequel we propose a projection onto these low- 
energy sectors which is artificial in a similar way as ND's 
splitting (||) in the one-band case but allows us to capture 
the important particle-hole excitations found within the 
perturbative approach. 



V. MODIFICATION OF THE 
NOZIERES-DE DOMINICIS APPROACH 

A. Splitting of the potential. Asymptotic forms of 
propagators 

Let us reconsider the academic problem above in which 
only the potentials V kk 2 — [V fc 2 ,y* are non-zero. To re- 
strict to energies close to /ii we assume these potentials 
to be separable as 



VJ 
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where the u k n are functions which refer to values of k 
in the subband n and which are maximum at the Fermi 
surface /i„/ and fall off sufficiently fast away from it. To 
maintain the hermiticity of the potentials V 12 we assume 
these functions to be real and equal, u]} = u 2 } = u\. 

In ND's approach, the constant part V 12 acts as an 
"external" potential on particles described by the opera- 
tors cjj = ^2 k u k c kn during times between t\ and t?,. The 
relevant physical quantities can be obtained from the so- 
called transient propagator, tp(t,t'), which describes the 
time evolution of these particles, with ti < t,t' < tg. 
The potential V 12 enters as the self-energy part into the 
Dyson equation for the transient propagator (see Eq. 
below) . The free part of this Dyson equation is described 
by the Green's functions 



k fe 

x [e{t)0{s kn - f x n )-0(-t)0{i,„ - .-,.„) . no) 



Following ND we must examine the asymptotic behavior 
of these functions as t — ► 00. We first consider the case 
n = 1. Here the cutoff function u\ is centered at the 
Fermi surface [i\ which corresponds to the situation of 
the original x-ray problem. Let us recall ND's arguments 
to obtain an approximate expression for G\{t): For large 
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times t the behavior of G\ (t) is dominated by the discon- 
tinuity at the Fermi surface. Let Qu{e) — i /1 (e)(u 1 (e)) 2 
with v 1 {e) the density of states in the subband 1, and 
^{ski) = u\. Then, for t > 



Gi{t) 



de gn(e)e 



(11) 



which becomes for large times t £ 1 (where £o is the 
characteristic energy describing the decay of u 1 ^)), 



Gi(t) 



-ilJ.it 



(12) 



The same expression holds for t < 0. We observe that 
the exact shape of u (s) is of no importance as long as 
\t\ ^> To include the short-time behavior, t < t^^ 1 , 
ND imposed that integrals over the product of G\ (t) and 
some slowly varying function (the transient propagator) 
must yield the correct result. Hence, if a is a time cut- 
off such that a ^q 1 but still much smaller than the 
characteristic time scale of the transient propagator, the 
short-time behavior of Gi(t) enters only through the in- 
tegral 



A x = 



/a poo | 
d«Gi(t)= / deg n (s)P 
-a J-oo £ - Ml 



= —Hi, 



(13) 

where we have used Eq. (tllj) for the second equality. This 
quantity can be added to the asymptotic expression jl^) 
in the form of A±5(t), and we obtain the central approx- 
imation of ND's approach 

G 1 (t) = -v n P±e^ t -9 n 8(t). (14) 

The case n — 2 is more subtle because there is no Fermi 
edge discontinuity at \i\ for electrons in the subband 2. 
Let us set g 2 \{e) = ^(^(u 1 (e)) 2 with v 2 (e) the density 
of states in the subband 2. We have, from Eq. (|i~0|), 



/>OC 

G 2 {t) = -i0(t) / deg 21 (e)e- 

J U9. 



(15) 



+ i8(-t) / deg 2 i(e)c~ 



(using Eq. (0)) 



-4, 



dtG 2 {t) =i dt deg 21 (e)e 



-%t(e— 



(17) 

The exchange of the two integrals demands some care be- 
cause the integrand becomes singular. For nonzero val- 
ues of e the time integration can be performed first and 
yields an integrand proportional to 1/e. In order to han- 
dle e — > 0, we notice that only the advanced part of the 
function survives because of the constraint t < 0. This 
implies that the energies entering in the exponential of 
Eq. ( |l7| ) must be shifted by an infinitesimal imaginary 
amount, e — > e + irj, with 77 > 0. The time integration 
can then be performed first, and we obtain 



A 2 = - 



/OO 
ds 021(e)- 
-00 ' 



1 



/ii + irj ' 



(18) 



where we have neglected the term at the lower boundary 
of the time integral since it is due to the sharp artificial 
time-cutoff a. With the formula l/(x — irj) = P(l/x) + 
iir5(x) we finally find 



-v 2 i + iirv 2 i 



with 



,00 1 

v 21 = / deg 21 (e)P 

■/-00 e- 



Mi 



(19) 



(20) 



which coincides with v 2 i, given by Eq. (^) with a slightly 
modified definition because of the weight function g 2 i(e). 

The approximate form resulting from ND's argumen- 
tation is therefore given by 



G 2 (t) = [i-Kv 2 i-V 21 ] S(t). 



(21) 



There is no algebraic decay in 1 /t because the Fermi sur- 
face discontinuity for the subband n = 2 is of no relevance 
for processes of energies close to e = fj,i . 



B. Transient propagators and absorption rates 



If fi 2 — /ii is large with respect to the decay of the func- 
tion u x (e) the first integral can be neglected. Using the 
same argument, we can push the upper boundary of the 
remaining integral to infinity. Hence 



G 2 (t) = +i9(-t) dsis^iu 1 ^)) 2 



(16) 



For large times, \t\ ^s> l;^ 1 , the period of oscillation of 
the exponential is much shorter than the time of vari- 
ation of the u 1 (t) so that the integral averages to zero. 
Following the above argumentation, the short-time be- 
havior of G 2 (t) can be resumed into a term A 2 S(t) with 



With the approximate forms of the initial Green's func- 
tions, Eqs. dTIf ) and (HJ), we can follow ND's analysis step 
by step, and we focus here only on the new features in 
the derivation. 

The computation of physical quantities is based on the 
transient propagators <p nn which can be defined by the 
Dyson equations 

V nn '{t,t') = G n (t-t')S nn , 

-J2 I * dt" G n (t - t")V nn " cp n " n ' {t" ,t'). (22) 
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For the academic problem under discussion, the functions 
(p 11 and (f 21 are found to form a closed set of equations, 



conditions of validity must be checked carefully at the 
beginning. 



At" G!(t-t")V 12 V 21 {t",t') 



ip 21 (t,t') = - / dt"G 2 {t-t")V 21 ip n (t",t'). 

(23) 

Since G2{t) is proportional to a (^-function, we obtain a 
single integral equation for ip 11 , 



V 11 (t 1 t') = G l {t-t')- / At"G 1 {t-t")V xl ^\t",t l ), 
in which the effective interaction 



-Vn = (-V 12 )[iTTv 21 -P21K-V 21 ) 



(25) 



leads to the same coupling constant g as in the perturba- 
tive approach, given by Eq. (|). The remaining problem 
coincides exactly with ND's original x-ray problem with 
the interaction potential V replaced by the effective po- 
tential V\i. Therefore, the absorption rate becomes, from 
Eq. (66) of Ref. |, 



A(ft) ~ Im 



I f2 — LJo\ 



with 



6 5 2 
e = 2 5 , 



(26) 



(27) 



and S the complex "phase shift" at e = fii induced by 
the effective potential Vn, 



tan S 



1 + v lx V\\ 



(28) 



For small V 12 , one has 5/ir ps v±iVi%, and the absorption 
rate merges with the parquet graph result of CR. Fur- 
thermore, the imaginary part of 5 leads to a broadening 
of the sharp Fermi edge singularities which exists in the 
case of disconnected subbands. 



VI. CONCLUSION 



The key approximation is the restriction to low-energy 
fluctuations close to the Fermi surface(s). The separa- 
tion of momenta in the interaction potentials must occur 
in these energy sectors. In the single-band case of the 
original x-ray problem, this low-energy restriction enters 
naturally by the discontinuity of the Green's functions at 
the Fermi surface which controls the large time behavior 
of these functions. The particular form of the weight 
functions u n {e) which project onto the low-energy sector 
has only little importance, so that it even might be taken 
as constant over the whole band, as often assumed in the 
literature. 

We have noticed with CR that the latter assumption 
leads to unphysical results within ND's approach in the 
case of nonequilibrium systems. Here the form of the 
weight functions becomes important, and we have shown 
how the approximations of the initial Green's functions 
must be chosen to obtain physical results. In particular, 
the decay in 1/t which is characteristic for ND's asymp- 
totic approximation, appears only if a Fermi surface lies 
within the energy range in which the cutoff functions 
u n (e) are non-zero. 

With this taken into account, ND's method remains 
valid and provides us with quantitative results for the ex- 
ponents of the power-law divergences at the Fermi edges. 
Moreover, the absence of the 1/i-decay in the subband 
n = 2 reduces the complexity of the discussed multi- 
channel problem. Instead of systems of singular integral 
equations, we deal with a scalar (single-channel) problem 
in which the presence of the other channels is reflected by 
a mere renormalization of the interaction strength which, 
however, becomes a complex quantity. 
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In our opinion, there remains an inconsistency with Ng's ex- 
tensions, independently to the topic discussed here, which 
arises from a too specialized Ansatz to the solution of the 
problem. A detailed discussion of this discrepancy is pro- 



vided in a forthcoming publication. The conceptual diffi- 
culty addressed in the present paper, however, persists also 
in a modified Ansatz. 



O 

FIG. 1: Two nonequilibrium Fermi seas extending between 
the energies (ei,/j,i) and (22,^2), respectively. A zero energy 
particle-hole pair between the subbands is indicated. Pairs 
with positive or negative energies are possible. 



